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An equivariant BRST-construction is used to define the continuum SU (3) gauge theory 

on a finite torus. I corroborate previous results using renormalization group techniques by 

explicitly computing the measure on the moduli-space of the model with 3 quark flavors 

to two loops. I find that the correction to the maximum of the one-loop effective action 

is indeed of order g 2 in the critical covariant gauge. The leading logarithmic corrections 

from higher loops are also shown to be suppressed by at least one order of g 2 . I therefore 

can relate the expectation value k 4 = — (Tr7 2 ) of the moduli 7 to the asymptotic scale 

parameter of the modified minimal subtraction scheme. An immediate consequence is 

the non-perturbative result (8 W ) = — gfs-ft 4 = — ^5/3^ Af^g = — 0-2228. . -A 4 -^ for the 

expectation value of the trace of the energy momentum tensor of QCD with three quark 

flavors. This relation compares favorably with phenomenological estimates of (O w ) from 
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QCD sumrules for the charmonium system and A4^= from r-decay. 
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1 Introduction 



The perturbative regime of ordinary SU(n) gauge theory with rif quarks in the fundamen- 
tal representation of the group has been thoroughly investigated in covariant gauges. The 
ultraviolet divergences of the loop expansion can be regularized by analytic continuation 
in the dimension of Euclidean space-time. The infrared divergences of ordinary covariant 
perturbation theory due to the massless gauge bosons and ghosts are however an indica- 
tion of the limited validity of this asymptotic expansion. The infrared problem is partially 
circumvented by the Operator Product Expansion (OPE) P] which is parametrized by non- 
perturbative matrix elements of certain local operators. Values for these matrix elements 
generally have to be extracted from experiment or from lattice simulations. I will here 
calculate one such non-perturbative quantity, the vacuum expectation value of the trace 
of the energy momentum tensor. The computation of this expectation value in terms of 
the asymptotic scale parameter requires little more than an asymptotic analysis of the 
theory. 

The basic idea for this calculation was presented in ref.i7|J. One uses the fact that 
infrared divergences are regulated in a gauge invariant fashion by considering the theory 
on a compact space-time manifold such as a torus. Asymptotic freedom asserts that 
perturbation theory is reasonably accurate for sufficiently small physical volume of the 
compact Euclidean manifold and one computes the deviations that arise as the volume of 
the torus becomes large [§. 

The main challenge is to define the model on a compact Euclidean space-time in a 
covariant manner. This separates (dynamical) degrees of freedom that propagate from 
those which do not (moduli). In a reasonably gauge- fixed theory on a torus, the latter 
are a finite set of parameters and the loop expansion for the dynamical degrees of free- 
dom depends on their values. In covariant gauges, that is gauges which preserve all the 
isometries of the compact space-time, the moduli are generally space-time independent 
constants which are related to vacuum expectation values of the basic fields and their 
composites. They differ from (dimensionful) couplings in that they take values in a cer- 
tain (moduli) space with a measure which has to be calculated. The volume dependence 
of this measure essentially determines the transition to the thermodynamic limit of the 
theory: at sufficiently small space-time volumes the measure on the moduli space is almost 
flat but it constrains the relevant moduli-space to a submanifold in the thermodynamic 
limit that is characterized by a few constants. 

The problem of constructing a sensible covariantly gauge-fixed SU(n) theory with 
quarks in the fundamental representation in a compact Euclidean space-time was solved in 
ref./|| using an equivariant BRST-quantization that eliminates the constant zero-modes 
of the Faddeev-Popov ghosts. The argument for this construction can be summarized as 
follows. With fermions in the fundamental representation of the gauge group on a torus, 
the gauge field and Faddeev-Popov ghosts in covariant gauges are periodic fields|]5|. 
Global gauge invariance of the action implies that the constant (anti-) ghost in ordinary 
covariant gauges decouples. One can show that this leads to a vanishing partition function 
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of such a continuum model. [One verifies]^, [|] that the Witten-index vanishes when the 
gauge fixing is viewed as a Topological Quantum Field Theory (TQFT) on the gauge 
group]. The equivariant BRST-quantization eliminates this problem without destroying 
the covariance of the gauge-fixed model. Generic zero-modes of the ghosts in this case 
are absent for the finite torus. The construction introduces constant ghosts with ghost 
number 0, ±1 and ±2 that form part of the moduli space. There are no zero-modes of the 
antiperiodic fermions to contend with and the remaining zero-modes of the gauge field are 
bosonic and do not give a vanishing partition function. These zero modes are addititonal 
moduli, but as observed in||, zero-modes of the gauge field on a torus with fermions in 
the fundamental representation are innocuous and can be neglected for sufficiently large 
volumes. 

Let me also justify considering the theory on a symmetric torus instead of some other 
compact manifold. The predominant reason is that this manifold preserves translational 
invariance and that the thermodynamic limit of interest is thus probably approached 
smoothly and uniformly. A torus of finite volume breaks the rotational invariance of 
Euclidean space, but this symmetry is not associated with any scale and generally believed 
to be recovered in the thermodynamic limit. A torus with periodic boundary conditions 
for the gauge fields also admits anti-periodic fermions in the fundamental representation, 
i.e. quarks. In addition, the space of gauge orbits on any finite torus with periodic 
boundary conditions is connected and the gauge dependence of topological sectors Q is 
therefore absent. Because the orbit space is connected, there is no strong CP- violation on 
any finite torus with periodic boundary conditions for the gauge fields. Instantons and 
anti-instantons appear (and disappear) pairwise and the Pontryagin number vanishes. 
Due to translational invariance there is some hope that the net effect of such meta- 
stable configurations can be summarily described by constant moduli. We will see that 
the moduli in fact saturate the vacuum energy density completely and that there are 
no "other" non-perturbative contributions to this quantity. A resolution of the [/^(l)- 
problem in this context is on the other hand still outstanding. [The solution proposed 
by Veneziano |7j appears promising if the axial vector ghost has an effective mass that is 
inversely related to the volume of the torus. The susceptibility would vanish for any finite 
torus but the {7a(1) problem could nevertheless be solved in the thermodynamic limit 
where this (conjectured) ghost becomes massless. Such an axial (Goldstone) vector ghost 
has however not been identified.] Let me note here that the {7a (1) "problem" is peculiar to 
covariant gauges, since its formulation requires the use of Goldstone's theorem to saturate 
Ward identities of the conserved, but gauge dependent, anomalous {7^(1) current. 

A technical argument for using a symmetric torus as compact Euclidean space-time 
is that perturbation theory on this manifold is relatively straightforward and that the 
thermodynamic limit of the expressions is readily taken. The mode expansion can be 
explicitly performed and an analytic continuation in the dimension D of the torus regulates 
the UV-divergent mode-sums of the loop expansion as it regulates the loop-integrals in 
ordinary (covariant) perturbation theory. 

In II the model proposed in I was qualitatively analyzed in detail. It was argued 
that perturbation theory for gauge dependent quantities is reliable in the vicinity of fixed 
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points of the parameter and moduli space only. Technically one observes that "dangerous" 
leading log corrections which tend to invalidate the perturbative expansion are absent at 
certain points of the parameter and moduli space. As we will see explicitly, the problem in 
perturbative QCD generally is that some of the leading log corrections of the loop expan- 
sion arise because generic gauge parameters and moduli flow to certain nontrivial fixed 
points. Using the renormalization group (RG) equation, these corrections can in principle 
be resummed (in practice this is approximately possible only if the gauge parameters and 
moduli are already sufficiently close to a fixed point), leading to a nonanalytic dependence 
on the coupling of the gauge parameters and moduli. The mess can be partially avoided 
by perturbatively expanding at the fixed points for these parameters from the outset. The 
remaining non-analytic dependence of the perturbative expansion is then due to the scale 
dependence of the loop expansion parameter g, asymptotically described by A asp, the 
only physical parameter of the theory. 

There is another (non-perturbative) argument for expanding at the fixed points of 
gauge parameters and moduli. For generic values of the parameters and moduli, the 
loop expansion is sensitive to small changes in their values. One may imagine that such a 
change could effectively be induced by taking certain non-perturbative configurations into 
account. An example in point is the scale dependence of the covariant gauge parameter a. 
Non-perturbative arguments [||] suggest that a is scale independent, since it is a coupling 
parameter of a TQFT. In a perturbative calculation, a does, however, generically show an 
asymptotic scale dependence. It is suggestive to consider this apparent scale dependence 
as due to Gribov copies of (perturbative) configurations that have been neglected in the 
loop expansion. At the fixed point, the combined effect of the neglected Gribov copies is 
at least asymptotically not scale dependent. Since this scale dependence of a is absent for 
any Green function, one may argue that the influence of other Gribov regions cancels at 
the fixed point and that the perturbative answer in this case is asymptotically correct. A 
similar argument can be used to infer that contributions from neglected non-perturbative 
configurations do not change the value of the moduli at fixed points. A perturbative 
evaluation of the model is thus asymptotically self- consistent only in the vicinity of a 
fixed point in the parameter and moduli space. The one-loop calculations in II revealed 
nontrivial fixed points in the moduli and parameter space of that model for nj < n, 
and the corresponding gauges were called Critical Covariant Gauges (CCG). The 1-loop 
effective potential for the moduli implied that the global SU (n) symmetry of the gauge- 
fixed theory is spontaneously broken to U(l) n ~ l in CCG and the (planar) contribution to 
the Wilson loop from the corresponding Goldstone bosons showed confining behavior. 

In section 2 of this article I give a concise description of the SU(n)-mode\ discussed 
in II and also present a possible extension for SU(n > 2). The action and moduli space 
one has to consider simplifies considerably in the special case nj = n. The CCG with 
5 = 1 , a = 3 for nj = n is a subset of the class of gauges with 5 = 1 in which some of the 
moduli decouple. The simplified action we will subsequently consider is given by (|2.13j ). In 
section 3 the two-loop effective action for the remaining moduli is calculated to two loops. 
It is expressed in terms of renormalized quantities in section 4 where I also argue that 
higher order corrections are suppressed in the CCG with a = 3. In section 5 the scale k 
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of the spontaneously broken SU (n) symmetry is finally related to Ajjg and section 6 uses 
this result to obtain the expectation value of the trace of the energy momentum tensor 
in terms of Ajjg. I conclude by comparing with the phenomenology of QCD sumrules. 



2 The Model 



I will consider a special case of the covariant model of Euclidean QCD on a finite D- 
dimensional symmetric hypertorus = L x L x . . . x L proposed in //. The boundary 
conditions for the gauge field A^x) and Faddeev- Popov ghosts c(x) and c(x) are periodic 
and the model is descibed by the tree-level action 



S = S c + 2 dxTr 



— (d-A(x)) 2 + Sc(x) D A -dc(x) + (1 - 5)c(x)d-D A c(x)+ 
2a 

1 , 1 



+a8(l — S)g 2 c(x)c(x)c(x)c(x) — S[c(x), c(x)Yf + - — -7 2 'jd-A 

2ag ag 

—ac(x)c(x) — ail — 5)<fic(x)c(x) -a<p H — jcfx) H — ac(x) 

9 2 9 9 



(2.1) 



where Sc is the Yang- Mills action of an SU(n) gauge theory with rif quark flavors in the 
fundamental representation 

i n f 
S c = [ -TiF^(A)F^(A)+J2^0- mi )^. (2.2) 

The Euclidean Dirac operator here is 2 

$^ i = l^ + gA^ l (2.3) 

and the field strength F^ U (A) is 

F^{A) = d„A v - d v A„ + g[A„ A v ] (2.4) 

Note that a CP- violating term proportional to the Pontryagin number j T Tr FAF vanishes 
on a torus with periodic boundary conditions^ I therefore did not include such a term 
in Q2.2|) . As discussed in the introduction, there is no strong CP violation on a finite torus 
with periodic gauge fields. 



The action is on-shell invariant under the equivariant BRST-symmetry generated 
by the operator s. Its action on the fields is 

sA^{x) = D A c(x)-[u,A^x)} 

2 The Euclidean Dirac matrices 7^ satisfy 7^7^ + 7^7/j = 2<5 M „1. I will generally suppress color indices 
for notational clarity. Except for ipi,Tpi, vectors of the fundamental representations of SU (n), all fields are 
traceless, anti-hermitian n x n matrices in this notation and transform under the adjoint representation 
of the group. The commutator [■, ■] is graded by the ghost number. 

3 "Toron" sectors with non- vanishing Pontryagin number exist on a torus with twisted boundary condi- 
tions. These however can only be imposed in the absence of quarks in the fundamental representation^]. 
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where 



sc (x) = -[u,c{x)\ - -[c{x),c{x)\ 

su = —~[w,u] + <f) 

s<p = —[u,(j>] 

sc(x) = —[uj,c(x)]+b(x) 

sb(x) = — [u, b(x)] + [(f), c(x)) 

scr = — [uj, cr] + a 

sa = — [uj, cr] + [(f), cr] 

S7 = — [u;,7] + 7 

si = -[w, 7] + [(f>,l] 

si[)i(x) = -uif)i(x) - gc{x)i)i(x) 

sijj^x) = -i)i(x)u - g4>i{x)c{x) 



1 

— ( 
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D*c(x) 



(2.5) 
(2.6) 



■g[A^(x),c(x)) 

is the usual covariant derivative in the adjoint representation. In (ZTT) the Nakanishi 
Lautrup field b(x) has been eliminated using the equation of motion 



b(x) = (j/g - d-A(x))/a - Sg[c(x),c(x)]. 



(2.7) 



It is straightforward to show that the BRST-operator defined above is nilpotent on 
any element of the graded algebra constructed from the fields of Table 1: 

s 2 = (2.8) 



field 


M x ) 


ipi(x)kipi(x) 


c(x) 


c(x) 


b(x) 





UJ 


a 


a 


7 


7 


dim 


1 


3/2 





2 


2 








4 


4 


2 


2 


(f)U 








1 


-1 





2 


1 


-2 


-1 





1 



Table 1. Dimensions and ghost numbers of the fields. 



The constant ghost u generates global gauge transformations of all the fields in the 
BRST algebra, except itself. The action does not depend on it and one restricts observ- 
ables to the equivariant cohomology S, 

dO 

S = {C: — = Q;sO = 0,O^sF} (2.9) 

where T is itself cj-independent. Since one is eventually interested only in expectation 
values of gauge invariant functionals of A, $ and ip, the notion of physical observables can 
be further sharpened to functionals in the equivariant cohomology with vanishing ghost 
number [|2], |!|. 
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For the sake of completeness let me note that ( |2.1| ) is the most general covariant 
power counting renormalizable classical action only for an SU(2) group. It depends on 
only two gauge parameters a and 5. For SU(n > 2) one can still extend the action by a 
BRST-exact and power counting renormalizable term proportional to 

sTr cA^A^x) = Tr b(x)A^x)A^x) - c(x)A^x)D*c(x) - A^{x)c{x)D^c{x) (2.10) 

The equation of motion of the Nakanishi-Lautrup field b(x) is then modified to 

b(x) = (j/g - d-A(x))/a - 5g[c(x), c(x)} - pgA 2 (x) . (2.11) 

Upon elimination of b(x) using fl2.11| ) the most general power counting renormalizable 
tree-level action of a covariantly gauge fixed SU(n > 2) model depends on three gauge 
parameters a, 5 and p and has the form 

Si = S + 2p f dxTig{d-A{x))A 2 (x)+g 2 a(5-l)[c(x),c(x)]A 2 {x) 

+gac(x){A^(x), ^c(x)} + ^pg 2 aA 2 (x)A 2 (x) - ^A 2 (x) , (2.12) 

where {•, •} denotes the anticommutator. For SU(n > 2) the analysis of II in this sense 
is incomplete and should be repeated with the more general action ( |2.12j ). Let us only 
note here that a non-trivial value of the bosonic ghost 7 would simultaneously regulate 
the infrared behavior of the tree level ghost- and gluon- propagators if a stable fixed point 
with 5^0 and p 7^ exists for the SU {n > 2) model. The additional parameter p could 
also soften the constraints found in 77 and lead to non-trivial fixed points of the moduli- 
and parameter- space of an SU{n > 2) gauge theory also for rif > n. Since 7 also couples 
to A 2 in (|2.12|) , the analysis of II becomes quite involved in this more general setting. 



My objective here is the relation between the perturbative asymptotic scale A^g^ and 
the fixed points of the moduli space in the simplest possible scenario and I will not 
embark on an analysis of the rather complicated model described by ( |2.12j ). For the 
purist the following considerations apply only to an SU{2) gauge theory, since the tree 
level action (ETTJ) in this case is the most general one. 



I will in fact restrict myself to the RG-stable subset of covariant gauges with 5 = 
l,p = 0. Using the result of //, one can expect a nontrivial fixed point in the moduli 
space of this restricted class of models for rif = n in a = 3 gauge. Choosing 5 = 1, p = 
greatly facilitates the calculations, because the moduli <p and a can be integrated out and 
there is no 7A 2 -vertex in this class of gauges. At 5 = 1 the coupling of <fi to [c(x),c(x)] 
in ( |2.1| ) vanishes, and the integration over the 0-moduli just sets a to zero. The quartic 
ghost coupling also vanishes at 5 = 1, further simplifying the calculation. 



For 5 = l,p = the tree- level action (2A) thus effectively is 

1 



S = S c + 2 f dxTr 



--(d-A(x)) 2 + c{x) D A -dc{x) - [c{x), c(x)] 7 + 7-^7' 
2a 2ag z 



•yd-A(x) H — 7c(x) H — ac(x) 

ag g g 



(2.13) 
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The perturbative expansion depends only on the moduli 7 of vanishing ghost number, since 
the moduli a and 7 just serve to eliminate the (on a finite torus normalizable) constant 
modes of the ghost and anti-ghost. The action fl2.13j ) is stable under renormalization and 
is the starting point of this calculation. 



3 The Measure on the Moduli Space 

The expectation value of an observable O G £ is formally given by the path integral 

(0)=Af J d^fdad-f J J [V^V^VAVcVc] O expS, (3.1) 



which defines the perturbative loop expansion. It was shown in / and II that ( |3.1| ) is 
generally normalizable, i.e. that (1) 7^ for an SU(2) group. The global bosonic ghosts 
7 introduced by the equivariant BRST-algebra are the only relevant moduli parameters 
of the model 0. We indicated in (|3.1|) that the integration over this finite dimensional 
moduli-space should be performed after the path integral over dynamical fields. As far as 
the dynamical fields are concerned, the moduli 7 are thus parameters of the action. On 
a torus, the integrations over constant ghosts a and 7 can be performed explicitly. They 
eliminate the zero-momentum modes of the Faddeev-Popov ghosts. The bosonic ghost 7 
in principle can also be eliminated, since S depends on 7 only quadratically. This would 
result in a nonlocal 4-point interaction of the dynamical ghosts that gives rise to infrared 
divergences in the loop expansion. A resummation of the perturbation series is therefore 
necessary and the term Tr 7[c(x), c(x)} in ( |2.13| ) should be treated as an unconventional 



mass term for the dynamical ghosts. The ghost "masses" , 7, are finally to be integrated 
with a certain measure. The measure, <i/i(7), on the moduli space 



n 

dfitf) = e rw J] dj a (3.2) 
0=1 

is described by the effective action T(7). We will calculate T(7) order by order in the 
loop-expansion for the dynamical fields. Note that T(7) is proportional to the volume of 
Euclidean space-time to leading order in L. In the thermodynamic limit the measure (|3.2| ) 
therefore effectively constrains the moduli space to the absolute maxima 7 of T(7)- The 
following calculation will determine this space. 



3.1 The Tree-level Measure 



To lowest order in h the measure on the moduli-space is Gaussian and apparently given 
by the quadratic term ^p-Tr 7 2 of the action ( |2.13| ). Note that the term / Tr7<9-yl in (|2.13| ) 
vanishes on a torus with periodic boundary conditions for the gauge field A^x). Since 
we eventually wish to consider the thermodynamic limit of the symmetric torus, 1/L is 
not a very useful mass scale. Introducing an arbitrary, but in the limit L — ► 00 finite, 
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renormalization mass /x, the tree-level effective action in D = 4 — 2e dimensions can be 
expressed in terms of the dimensionless quantities that also appear in the 1- and 2-loop 
contributions 



where g = gfi~ e is the dimensionless coupling constant and the e,,i = 1, . . . ,n, are the 
real eigenvalues of the traceless hermitian matrix i'j = i^ a t a . Since Ya e % = 0? the effective 
action is a function of the dimensionless differences 



% - ^ • (3-4) 



of the eigenvalues only. 



Note that (|3.3|) is the effective action in terms of bare parameters a = Z a aR, g = Z g g R 
and 7 = Z^ R . By calculating the (regularized) effective action to two loops, I will 
determine the first few coefficients in the loop-expansion of some of these renormalization 
constants in the MS'-scheme. The coupling g and gauge parameter a are also related to 
vertices of the dynamical fields and one can thus check the validity of Ward-identities to 
some extent explicitly. 

Before calculating loop corrections to the effective action note the following interesting 
fact. Let us for the moment ignore the constraints imposed by the BRST-algebra ( |2.5| ) and 
regard 7 as the zero-momentum mode of a background field 7(2;). The coupling <y{x)d ■ 
A(x) of j(x) to the longitudinal gauge field in (|2.13|) in this case precisely compensates 



the 7 2 (^) _ t erm °f the tree- level effective action for all non-zero momentum modes of 
j(x). The effective action for non-constant modes of 7(2;) in fact would vanish to all 
orders in the loop expansion due to a supersymmetric compensation between ghosts and 
longitudinal gluons. This supersymmetric compensation is absent for the constant part 
7 of j(x). The effective action for constant 7 does not vanish, and the corresponding 
measure on the moduli space is therefore not flat. In this sense the non-trivial measure on 
the moduli space we will find is a consequence of the lack of supersymmetric compensation 
between ghost and (unphysical) gluonic degrees of freedom with vanishing momentum. 
The same imbalance ensures that the partition function, proportional to the Witten index 
of the TQFT on the gauge group p, |J, does not vanish. This mechanism is reminiscent of 
the lack of supersymmetric compensation for the ground state in models with unbroken 
supersymmetry and thus perhaps worth noting. The following will determine the ground 
state configuration of the interacting theory. 



3.2 The One-Loop Measure 

The measure on the moduli space is Gaussian to lowest order in H only. The dynamics 
depends on the value of 7 and the one-loop radiative corrections to the effective action 
are depicted in Fig. 1. These contributions to the effective action were evaluated in iT 
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(and for SU(2) also in0). For a D = 4 — 2e dimensional torus, the leading contribution 
for Ljj, ~ oo is0 

r^(7) = (»L) D £ 2«S-oos(i7re)r(e-2). (3.5) 

l<i<j'<n 



2p 



7 r 



Legend : 



- "V ..J._ y- insertion 

ghost-propagation 

Fig. 1: Diagrammatic representation of 1-loop contributions to the effective action of the moduli. 



Defining 



tl3 := lE -l + Un(vl) , (3.6) 



where 7# = 0.577216 ... is Euler's constant, fl3.5|) for e ~ has the expansion 



l<i<j<n 



Note that r 1_ ' oop does not depend on a nor on g. Since 7 = Z^Tb with — 1 of O^g^/e), 
pi-ioop gj ves r j se ^ a t erm f order ^ in the renormalized two-loop effective action. This 
contribution is proportional to 



= W D E (-2)«.) 1 - s/2 cos(^vr)r(e-l) 

8=1 l<«<i<n 

= W D e < f--2^ + e[l-^ + 4] + 0(, 2 )) (3. 

l<i<i<n iZ / 



where the term of 0(e) still gives a finite contribution to the renormalized two- loop 
effective action, since the coefficient — 1 is of 0(g 2 R /e). 



3.3 The Two-Loop Measure 

The 7-dependent part of the (unrenormalized) 2-loop effective action is depicted in Fig. 2. 
With the methods of //, the leading contribution for [iL ~ 00 can be computed for a 
/^-dimensional torus. The zero-modes of the ghosts are absent in the mode expansion 
and the zero-modes of the gluon field give a subleading contribution for large fiL. For the 
leading contribution to the effective action the mode-sums can be replaced by integrals 
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since they are IR-finite in D > 2 for 7 7^ 0. In appendix A these integrals are evaluated 
using conventional Feynman techniques with the result ( |8.11|) , valid for 2 < D < 4. For 
e ~ the expansion of ( |S.11|) has the form 



r 2-loo P( - A e~0 ( n D n<3 2 V- 2 f 3-Q 2 - (3 - a) [1/4 + tjj] 

r p ( 7 ) = (/iL) — - ^ u . \^- 2 ~+ + 3-4ty 

+ (a - 3){finite terms} + 0(e) \ . (3.9) 



V v / 



q 



r" 



*AAA/\A/WW\> = SAAAAAAAAAAfc - ^AAAAAAAAAA/* 

"1 p-q £- *. p-q # \ p-q / 



n-l x >. n .-C Sj, P y \ n 

• • • 

2n~ 

a) b) c) 

Fig. 2: a) Diagrammatic representation of the 7-dependent 2-loop contributions to the effective 
action of the moduli, b) The contribution with the tree-level ghost propagator (hatched) for 
fixed moduli 7^0. c) The 7- independent subtraction at 7 = 0. 

Note that the terms of ( |3.9| ) proportional to %/e and 1/e 2 vanish for a = 3. The 
two-loop contribution to the effective action is therefore of order g 2 instead of order 
gHij compared to the one-loop effective action ( |3.7|) in this gauge. We will see that this 
suppression of the two loop corrections in a = 3 gauge is not accidental and that higher 
order terms of the form ^ 2n s™ +1 with n > 1 are also absent in this gauge. 



4 The Renormalized Two-Loop Effective Action 

We have so far computed the dependence of the effective action on the bare variables 7, g 
and a in D < 4 dimensions to two loops. To take the limit e — > 0, the effective action 
has to be reexpressed in terms of renormalized quantities that remain finite in this limit. 
We could of course have computed T as a function of 7r> 9r an d otR from the outset, 
by writing the action (|2.13| ) in terms of renormalized fields and couplings and including 
counterterm contributions in the computation of the effective action. The two methods 
are entirely equivalent, but the approach I chose requires less computational effort. It 
avoids the computation of an infinite set of (counterterm) diagrams which just give the 
derivative ( |3.8| ) of the infinite set of 1-loop diagrams already computed. 

To express the two-loop effective action by the renormalized quantities, we simply 
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substitute the expansions 



7 


= IrZ^ 


= Jr( 


1 + + 


f 


- a 2 7 2 

- 9 R ^g 


= 9l( 


1+ & z« + 

16n 2 


a 


= OL R Z a 


= an 


{ + 16vr 2 a 



Z? ] +0(g R )) 



(167T 2 ) 2 

-^V^ 2) +0(^)), (4.1) 



(16tt 

for the bare quantities in in the previous expressions. In the MS'-scheme the are 
determined by requiring that they cancel only the l/e k terms. To the order we calculated, 

T 2-Ioo P{ - rZ _. *a aRZa) = T ^oo P{ - r . ^ aR) + q ^ r) ^ (4 2) 

and 

+ 0(&). (4.3) 



=i 



Comparing ( |3.9| ) and fl3.8| ), the vfjtij/e term in ( |3.9| ) requires that we choose 



Z M = 3 ) (4.4) 



in the MS'-scheme. The remaining divergent terms of, 



^ + ^ = {flLf ^ Jaj (I + M (<!_-) + <L£*) ) + flmte tems 

(4.5) 

have to be canceled by the tree-level effective action ( |3.3| ) expressed in terms of j R , a R 
and g R 

if? 2 2 / -2 
l<i<i<n na R9 R \ iD7r 

+ (J^V ((4 1} - - Z^f - Z (1) Z ( ^ + 2Zf - Z^ - Z^)\ + 0(g R ) (4.6) 



Comparing with ( f4.5| ) one obtains two relations for the Z's in ( |44] 



z (i) + z (i) = 2 ^ 1} - — = n<y0lR + ^ 



2e 

2Zf = Z (2) I Z {2) I Z (1) Z ( a 1] n2 (^ a R + ^R + 9 ) , n 2 a R (5 + a R ) 

The first of these is a check of the two-loop calculation, since Z^ l > and Z^' can also 
be determined from the 1-loop radiative corrections to dynamical vertices. The second 
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relation in (4~7?) gives in terms of the renormalization constants for the couplings. 
With[| 

Z (l) = -fl,/e = (|n, - yn)/e 

= i^^n-\n f )/e (4.8) 

one verifies that the first relation in 04.71) indeed holds. 



The terms proportional to 1/e and 1/e 2 thus are seen to vanish through order g\ when 
the effective action is expressed in terms of 7^, §r and ocr. In the MS'-scheme we have 

T(7i?., 9r, oir) = (^L) D v Rij { — + V 2 - ^ffii 

l<i<j<n I na R9R 

+tH ( 3 / 2 " 2t »i + (<*R ~ 3){finite terms})) + 0(g R ) (4.9) 

l07T z J 

To simplify expressions allow me to sometimes abuse notation in the following and denote 
renormalized quantities without an index (R) whenever confusion with bare quantities is 
unlikely. 



4.1 Higher Loop Corrections to the Effective Action 

Since the coefficient given by ( J4.4Q , vanishes for a R = 3, the anomalous dimension of 
7r is of order in this gauge. In II it was observed that = 3 is the stable non-trivial 
UV fixed point of the gauge parameter for nf = n because given by (|4.8|) , then 
vanishes as well. Using RG-arguments, it was shown|2| that a = 3, 5 = 1, (p = 0) is the 
CCG for rif = n. Remarkably the leading asymptotic correction of order h in the effective 
action ( |4.9|) is of the form g^tmj rather than gRS 2 Ri j for ocr = 3. The two-loop contribution 
to the effective action in this gauge is thus suppressed by a factor g 2 compared to the 
one-loop contribution. We will see explicitly that this is no accident and that in fact all 
higher loop contributions are similarly suppressed asymptotically in a R = 3 gauge. 

It is not difficult to show that a term of order 7^^ in the renormalized effective 
action is associated with an unrenormalized fc-loop contribution that diverges as 7 2 /e k . I 
will now argue that a generic 1-Particle Irreducible (1PI) k > 2-loop vacuum diagram that 
depends on 7 at most diverges like 7 2 /e fc_1 for a = 3. The vacuum diagram we consider 
obviously must have at least one ghost loop and there must be at least two insertions of 
7 for it to depend on 7. We select a ghost loop with at least one insertion of 7. For n > 2 
the diagram is then generically of the form shown in Fig. 3, since a gluon must be emitted 
somewhere before the insertion of 7 on the ghost loop and eventually must be reabsorbed 
by the same ghost loop if the diagram is 1PI. 
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Fig. 3: Diagrammatic representation of the sceleton of a 1PI contribution to the 7-dependent 
effective action with two or more loops. Shaded areas denote generic parts of a multi-loop 
diagram that can be considered as a correction to the ghost-gluon vertex or of the two-point 
functions themselves. See the text for details. 

The leading divergence of all the k—2 loop integrals apart from li and I2, is furthermore 
at most l/e k ~ 2 . The main point of the argument is that such a divergence (if it occurs) 
cannot depend on 7 and is compensated by counterterms for the subdiagrams that do 
not depend on the mass-scale 7. The reason is that an insertion of 7 renders the ghost- 
gluon vertex and ghost selfenergy superficially convergent (since one of the ghost momenta 
factorizes in the gauge 5 = 1). The insertion of 7 thus reduces the leading divergence of 
these sub-diagrams by at least one power of e. A similar argument holds for the gluon 
polarization, since a single insertion of ('j)ab = f abC/ y c gives a vanishing contribution due 
to Bose symmetry. As far as the leading divergence of order l/e k ~ 2 for the subdiagrams 
is concerned, the diagram of Fig. 3 is thus equivalent to the two-loop diagram of Fig. 2 
multiplied by (const. /e k ~ 2 ). We have explicitly shown that the two-loop contribution 
( |3.9|) diverges only as 7 2 /e when a = 3. The leading divergence of the k > 2-loop 
contribution is thus at most 7 2 /£ fe_1 . The leading dependence on ln(// 2 ), viz. tjuj, of the 
renormalized effective action in D = 4 dimensions evaluated to iV-loops in an = 3,5 = 1 
gauge in the MS'-scheme is therefore of the form 

rW(,) = £ £J ei ~ e>? B£ + 5- l " + t Ck • (410) 

where the Ck are /^-independent coefficients determined by the (k + l)-loop contribution, 
ei,i = l,...,n are the n real eigenvalues of the traceless hermitian matrix and 
9 — 9ms is the renormalized coupling constant in the minimal subtraction (MS) scheme. 
For g 2 tij of 0(1) asymptotically, contributions to the effective action from two and more 
loops are thus suppressed relative to the one-loop contribution by at least one order in 
g 2 . The perturbative calculation of the effective action in the limit g — > thus gives a 
sensible asymptotic expansion in the gauge a = 3, 5 = 1, p = 0. 
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5 The RG- invariant maximum of the effective action 



Since the effective action r(7) in ( |3.2|) is essentially proportional to the space-time volume 
for (/uL) 4 ~ oo, only absolute maxima 7 of IX7) are of relevance in the evaluation of 
expectation values ( |3.1| ) of observables in the thermodynamic limit. We now determine 
the absolute maxima 7 of T(7) in a = 3, 5 = 1 gauge. 

In general the eigenvalues of 7 will be functions of the renormalization point /i, and 
will either vanish or approach infinity as — > 00. For nj = n, they however approach 
a finite value proportional to Ajjs- Maxima 7 of the effective action (|4.10| ) for large 
\j? I (e, — ej) satisfy 

8vr 2 d 



=1 



L 4 ds 

n2 J 16vr 2 t cig 

l<j<j'<n 



2 . „2 iV-2 



fc=0 



(5.1) 



where 



c fc (/u) = c fc+ i + -g 2 (k + 2)c fc+2 (5.2) 

and g 2 is the running coupling constant of the MS-scheme. Its dependence on the renor- 
malization point ji to two loops for sufficiently large fi/A MS is|| 

9 2 1 A lnln/i 2 /A|, /5 



16tt 2 A) ln/i 2 /AL s A) (A) ln/U 2 /A 2 //s ) 2 ' 
where 

lln — 2n/ , . 34n 2 - 13n f + 3n//n 
A) = 3 and /3i = ^ ^- (5.4) 

are the first two coefficients of the /^-function and Ams is the so defined /i-independent 
asymptotic scale parameter of the MS'-scheme. For rif < yn quark flavors g 2 vanishes 
logarithmically as /1 — > 00, implying asymptotic freedom. For fl5.1|) to hold in the limit 
/Lt — > 00, tjj must be of order 1/g 2 and we see that the sum in ( |5.1| ) is a correction of order 
g 2 to ty. Since the perturbation series diverges for N — > 00, this is more precisely true in 
the sense of an asymptotic expansion: for any fixed, but not necessarily small, order in 
the loop expansion, the contributions from two and more loops become negligible in the 
vicinity of the maximum ~ of the effective action as \x — ► 00. 



The condition ( |5.1|) for the maximum of the effective action can be regarded as an 
asymptotic definition of the coupling g in terms of the renormalization point fi and the 
scale k 4 = — Tr7 . The situation is similar to comparing the asymptotic scales of two 
perturbative renormalization schemes 0. Using ( |5.3|) to eliminate the coupling g 2 in favor 
of the asymptotic scale Ams of the MS'-scheme, ( |5.1|) in the limit \i — > 00 implies that 




... = (3 \n^/A 2 MS , (5.5) 

2^ ( e i ~ e i 



l<i<j<n 



14 



Eq. ( p.5| ) gives the asymptotic scale dependence of the eigenvalues e"j(//) corresponding to 
a maximum of the effective action. We see that the maximum approaches a finite limit 
as /i — > oo only if 3n = /3q- From (|5.4j) , this is the case for 

rif = n (5-6) 

quark flavors. 

For an SU (2) group and two quark flavors, (|5.5|) can be explicitly solved for the one 
independent eigenvalue e = e i = — e 2 of 17, 

e = ±2ne l -^A\ IS = ±|Aj^ ~ ±(1.1658A I7 ^) 2 (5.7) 

We thus have that 

k\u = nj = 2) = -Trf = 5? + = yA 4 ^ ~ (1.3864A^) 4 (5.8) 



In the case of SU(3) with three quark flavors, (|5.5|) by itself is not sufficient to 
determine the texture of the eigenvalues of 17. However, if the maxima of ( 4.1(jp are 



unique up to permutations of the eigenvalues e^, the discrete symmetry 7 — > —7 of the 
effective action must be equivalent to a permutation of the eigenvalues e"j. For an SU(3) 
group this requires that the absolute maxima of the effective action are permutations of 

ei = -e 2 = e, e 3 = . (5.9) 

For rif = 3, ( |5.5| ) then determines the value of e in terms of A]^ 

e = 2- 2 / 3 e(4vrAi f5 e-^) = 2^ 2 / 3 eA|^ ~ (l.SOSGA^) 2 . (5.10) 
and the corresponding value for the scale 

K 4 (n f = n = 3) = -Tr7 2 = 2- 1/3 e 2 A|^ ~ (1.5562A I7 ^) 4 (5.11) 



For SU(n > 3) groups and rif — n, symmetry arguments and (|5.5|) alone are not 
sufficient to completely determine the absolute maxima of the effective action. The re- 
lation ( |5.5|) however implies that one can find the fixed points in the limit \x — > 00 by 
(numerically) maximizing 

^^-Si^-^Tfe^- (5 - 12) 



This is not to say that ( 5.12| ) is the effective potential for 7. It cannot be, since the 
anomalous dimensions of 7 and a only vanish to order h for a = 3 and rif = n. T e ^ 
therefore does not satisfy the correct renormalization group equation. We have however 
seen that higher loop corrections become negligible in the vicinity of maxima for /j,/Aj^ ~ 
00, since Q5.5D shows that they are of order g 2 . Y e ^ correctly describes the variation of the 
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effective action near its maxima in the limit c? 2 — ^ and thus correctly reproduces 
in this limit. Higher loop corrections to ( [4.1 0|) in general are not negligible sufficiently 
far from the maxima when |ln |(e* — ej)//i 2 )|| ^> 1/g 2 . T e ^ of (|5.12|) in this case differs 
significantly from the true effective action. 

It was argued in II that the deviation of the true effective action from Y e ^ is of no im- 
portance in the thermodynamic limit, since the measure on the moduli space, proportional 
to e v ^\ constrains the moduli space of 7 to the immediate vicinity of absolute maxima of 
the effective action. Fluctuations of the moduli vanish as 1/L 2 . In the thermodynamic- 
(L — > 00) and critical- (g 2 — > 0) limit, the measure on the moduli space is thus effectively 
given by T e ^ . 

The following remark is perhaps of interest: our perturbative evaluation is uncontrolled 
near 7 = and the fact that r^(0) and r e -^(0) vanish does not imply that the true 
effective action also vanishes for 7 = 0, since |ln | (e< — ej)//i 2 || — > 00 1/g 2 for any 
finite value of /1 in this case. Higher orders in the loop expansion are therefore much 
more important than the computed ones and an asymptotic expansion of the effective 
action is pretty useless near 7 = 0. One could imagine that after a resummation of the 
perturbative expansion, the factor 7 2 of the effective action cancels for 7^0 and that 
the true effective action has a finite or perhaps even divergent limit for 7 — >• 0. In II it 
was argued that an absolute maximum of the true effective action at 7 = would imply 
that nonperturbative effects dominate even asymptotically, since the absolute maximum 
of the effective action sets the asymptotic physical scale of the model - as we have seen, 
this maximum does not occur at 7 = in any finite order of the loop expansion. This 
argument however does not a priori imply that the true effective action vanishes at 7 = 0, 
i.e. that r e -^(7) not only reproduces the true effective action near its absolute maximum, 
but also gives the difference in energy density between the broken and unbroken phases. 
The next section shows that r e " (7) does indeed reproduce the correct energy difference. 



6 (O/m) in terms of A 



We can also compute the energy density difference between the broken phase with (7 2 ) = 
7 2 7^ and the S'?7(n)-symmetric one with (7 2 ) = by a change of scale, since Euclidean 
invariance relates the vacuum expectation value of a component of the energy momentum 
tensor to the vacuum expectation value of its trace 

(e^> = 5 -f (e PP ) (6.1) 

and an insertion of W is the response to a change in scale. In the chiral limit with 
irtf = the action ( |2.13| ) in D = 4 dimensions is invariant under an infinitesimal dilation 



and a corresponding change of integration variables, 

Sx = —x5X =>5L = —L5X 
6$i(x) = (diQiW-Xud^xfiSX (6.2) 
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where $j is a generic field variable and di is its canonical dimension as given in Table 1. 
Note that the change of variables also includes the moduli, which are integration variables 
and in this sense differ from mass parameters and coupling constants. The reason for the 
formal invariance of the massless theory under a change of scale of course is that there are 
no dimensionful parameters in D = 4. This is no longer true for the regularized model 
in D < 4 (or for that matter any other regularization of the model) and ( |6.2| ) is not a 
symmetry of the renormalized theory (which does depend on a physical mass scale such as 
that does not vanish in the chiral limit). One can compute the scale anomaly from 
the fact that the (bare) coupling constants and fields of the dimensionally regularized 
model in D < 4 effectively change when the fields are transformed according to (|6.2|) . 
This anomalous contribution to the trace of the energy momentum tensor of QCD was 



first derived by Collins, Duncan and Joglekar|T(J. They found that the dilation (|6.2|) is 



equivalent to a zero-momentum insertion of the renormalized operator 

ev = + £ + im{g)) (6.3) 

in gauge invariant correlation functions. In ( |6.3| ) (3(g) is the /3-function describing the 
scale dependence of the coupling constant and 'jm(g) is the anomalous dimension of the 
quark masses. Naively one expects only an insertion of J2f m f^ f m the massive case. 
The rest of ( |6.3| ) is the anomalous contribution to the trace of the energy momentum 
tensor. The arguments of |10| also hold for the action ( p,13|) and the trace anomaly 



need not be recalculated here. It suffices to note that the action (|2.13|) differs from the 
model considered by[T(J in the gauge fixing sector only. Since BRST-exact insertions in 



gauge invariant correlation functions vanish in any dimension, the net effect of a dilation 



on gauge invariant observables is an insertion of ( |6.3j ) also for our model. Taking the 
vacuum expectation value of ( |6.3D 



= ^ (Wu) + £«V (*/*/) (i + iM) (6-4) 



we see that non- vanishing gluon and quark condensates imply a reduction of the vacuum 
energy compared to a hypothetical scale invariant phase in which these condensates vanish. 



Since a dilation (|6.2|) is equivalent to an insertion of MAl at zero momentum, we can 
compute ( |6.4| ) from the effective action IX7) for the constant ghost 7 after all dynamical 
fields have been integrated out. The finite dimensional integral over the moduli 7 itself 
does not diverge and the limit D — > 4 of the (renormalized) effective action can be 
taken before performing these integrals. The trace of the energy momentum tensor is 
furthermore RG-invariant and therefore does not depend on the renormalization scale \x. 
We can take the limit fi — > 00 to compute it. In this limit the change of scale (|6.2| ) 
effectively results in an insertion of — L A nTr 7 2 /(87r 2 ) since r 6 "^) describes the measure 
accurately in the vicinity of a maximum of the effective action. Since the anomalous 
dimension of 7 and a vanish to order g 2 in a = 3 gauge, we obtain for n/ = n flavors, 
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Ir'v = 

8vr 2 1 8vr 2 



Tr7 = — — -k (jif = n) (6.5) 



The expectation value of the trace of the energy momentum tensor can thus be ex- 
pressed in terms of the fixed point 7 of the maximum of the effective action r(7). Since 
the left hand side of (|6.5|) is the expectation value of a gauge invariant operator, it is 
worth emphasizing that the fixed point 7 of the effective action is not gauge dependent 
although our perturbative evaluation of 7 in terms of A;^ for rif — n flavors apparently 
was restricted to the gauge (Xr — 3, 6 — 1, p — 0. As emphasized in II, = 3 is however 
not just any particular gauge. It is the stable fixed point of the gauge parameter for 
rif = n. As shown in // this is related to the fact that Z a = 1 + 0(g A ) at «r = 3 for 
rif = n. In this sense lim^oo oir(ij) = 3 for whatever gauge (an 7^ 0) one starts from at 
finite [i. [To have lim^oo afi(/i) 7^ 3 for rif = n, the bare gauge parameter as(/i) would 
have to be /i-dependent - invalidating the RG-analysis. Note that Landau gauge, an = 0, 
is an unstable fixed point for rif < ^ n 0-] The "dangerous" leading logs in gauges 7^ 3 
are related to the running of o>r and the anomalous dimension of 7 in these gauges. In 
the limit — > 00 in which a_R(/i) — > 3, the contribution to the fixed point of the effective 
action from "dangerous logs" cancels (which can be inferred from the RG-equation) . In 
a loop expansion of the effective action to finite order this cancellation is not seen in 
gauges aji(fi) 7^ 3 because the perturbative expansion is analytic in g 2 . Choosing the 
fixed point an(fi) = 3 is not a choice of gauge in the fi — > 00 limit, but rather presents 
a simplification of the asymptotic analysis. The value a^(/i) at finite /i only determines 
how the fixed point an = 3 is approached [§]. 

The asymptotic evaluation of the expectation value of the trace of the energy momen- 
tum tensor in terms of k for nf = n flavors in ( |6.5| ) thus does not depend on the choice 
of covariant gauge at finite /i (with the possible exception of Landau gauge, an = 0, for 
which our asymptotic analysis is not valid). We can combine ( |6.5|) and ( |5.11|) to the 
physically interesting relation 



<e w > = -^\n f = n = 3) = -^Afe ~ - (0.687A^) 4 (6.6) 



_3e_ 2 

between the vacuum expectation value of the trace of the energy momentum tensor and 
the asymptotic scale parameter A-^=r of the modified minimal subtaction scheme of an 
SU(3) model with three quark flavors. 



Using ( |6.1|) the difference Ae between the energy density e true of the true ground state 
and that of the perturbative vacuum with unbroken S'f/(3)-symmetry is 

As = e tmc -e pcvt = i ((6 W ), - (6 W ) ?=0 ) = -7^^ = -yJ^Afe ~ -(- 486A ms) 4 - 

(6.7) 
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for SU(3) with three quark flavors. The difference Q6.7|) is precisely T e ^( / ~f) — T e -'-'(0) 

of (CT). 



7 Discussion and a comparison with phenomenology 



I wish to stress that the final result ( |6.6| ) is a determination rather than an estimate of 
the expectation value of the trace of the energy momentum tensor for an SU(3) gauge 
theory with three quark flavors. It implies that the unbroken phase with global SU(n) 
color symmetry, and thus 7 = (7) = 0, corresponds to the scale invariant phase with 



tyftyf) = and (F® V F^ V ) = 0. The relation (|6.q) furthermore shows that there are no 
other non-perturbative contribution to the expectation value of the trace of the energy 
momentum tensor. More precisely: if a nontrivial fixed point 7 of the moduli 7 exists, its 
magnitude gives the expectation value of the trace of the energy momentum tensor. As 
I have argued, this fixed point of the moduli space in fact does not depend on the gauge 
parameter a. Eq. (|6.6|) gives it an even broader, completely gauge invariant, meaning. The 
contribution to the vacuum energy density from non-perturbative field configurations such 
as instantons, monopoles and whatever else there might be is subsumed by the expectation 
value of the moduli 7- The perturbative expansion around a particular ground state with 
7^0 makes sense because Eq. ( |6.6| ) implies that this vacuum has the correct energy 
density. It is furthermore clear that the modified perturbative expansion automatically 
includes power corrections to physical correlation functions at large Euclidean momenta. 
These power corrections first appear at the three-loop level in correlation functions of 
gauge invariant quark currents. They arise due to insertions of 7 i n ghost loops. Due 
to global gauge invariance of physical correlation functions, the power corrections in the 
chiral limit are of order 7 /p 4 ,7 /p 6 , • • -etc. as one also expects from the OPE[jT|. The 



normalization of the coefficients is, however, no longer arbitrary. Eq. (|5.9|) and (|5.10| ) 
relate the asymptotic expectation values of 7 to the fundamental scale Aj^g of the model. 
We furthermore know that corrections to these asymptotic values are analytic in g 2 in 
the gauge a = 3, 8 = 1, p = 0. They can thus be reliably computed in the framework 
of perturbation theory. In this asymptotic sense the SU(3) model with nf = 3 flavors 
is solved, since the non-analytic "non-perturbative" corrections of the OPE have been 
determined in terms of the asymptotic scale parameter. One could consider this an 



explicit realization of the program first proposed by Stingl[ll 



We have here obtained the fixed point of the moduli space only for the special case 
of an S"?7(n)-symmetric model with nf = n flavors. The analysis in II suggests that a 
fixed point 7 7^ can also be found for nf < n. The perturbative analysis in this case 
is complicated by the fact that the CCG occurs for 5^1. The moduli and a then 
have to be taken into account and cannot be trivially integrated out. A non-trivial fixed 
point of the moduli space of an S77(3)-model with rif ~ 6 > 3 quark flavors, if it exists, 
would be of greater phenomenological interest. In the restricted space of covariant gauges 
with p = no such fixed point was ascertained^]. We are currently investigating whether 
there is a fixed point in the moduli space of the SU (3)-model for nf > 3 in the enlarged 
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set of covariant gauges with p ^ 0. 



By performing a two-loop calculation of the effective potential, we have explicitly 
verified that there is a non-trivial fixed point 7 in the moduli-space of an 5'[/(3)-model 
with nf = 3 quark flavors. The eigenvalues of 27 have been related (see ( |5.9|) and ( |5.10|) ) 
to the asymptotic scale parameter Aj^. The reader probably noticed that the one-loop 
effective action actually suffices for this determination. The two-loop calculation in this 
special case, however, explicitly demonstrates the suppression of higher order corrections 
to the fixed point of the moduli space in CCG - which inH was inferred from the RG- 
equation. The methods of appendix A can furthermore be used to calculate the loop 
expansion of S'f/(n)-invariant (physical) correlators in terms of the eigenvalues of the 
fixed point ij. Although these calculations are still somewhat cumbersome, we can in 
principle calculate all power corrections due to the moduli in this model order by order 
in the modified perturbative expansion. 

A comparison with phenomenology is complicated by the fact that more than three 
quark flavors are dynamical at momentum scales above /i 2 ~ m? c ~ (2GeV) 2 . An 577(3)- 
model with only three flavors can thus only be expected to have reasonable accuracy 
at energies where the effective coupling ot s is still rather large, a s (2GeV) ~ 0.35. Let 
us nevertheless attempt a comparison with the phenomenology of QCD sum-rules while 
keeping this restriction in mind. 

From moment sum rules in the J/\l/ channel and chiral perturbation theory, the LHS 
of ( |6.6D can be estimated |12|. Any reasonable |L2|, [13], [TJ] value for the gluon condensate, 
(tyty^ and the light quark masses gives that the quark contribution to (6 W ) is at most 
20%. It is thus comparable with the systematic error in the determination of the gluon 
condensate. For a rough comparison with QCD-sumrules (and that is unfortunately all 



we can hope for), ( |6.6| ) essentially states that 



(3) 

where A.jJ s is the asymptotic scale parameter of the modified minimal subtraction scheme 
for a three-flavor model. The best determination of the gluon condensate at the highest 
energies for which only three flavors are dynamical is obtained using QCD-sumrules for 
the charmonium system. These balance perturbative contributions against power correc- 
tions. Using ratios of moments, Shifman, Vainshtein and Zakharov (SVZ) [|T2]] originally 
estimated 

—KuKu) ~ (330 ± 30MeV) 4 with a s (2.5GeV) = 0.2 and m c ~ 1.26GeV 

7T M ^ / SVZ 

(7-2) 

Most recent high-energy data indicate that their value for a s is a bit low and that 
a s (2.5GeV) ~ 0.3 would be more appropriate (see Fig. 4). Since a s determines the 
perturbative contribution to the moments, the extracted condensate value depends some- 
what on a s . In the method used by|Lj] it is also relatively sensitive to m c . 
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Fig. 4: The running strong coupling constant a s of the modified minimal subtraction (MS 1 ) 
scheme. The experimental data and best fit (solid line) with its la error (grainy area) are 
reproduced from figure 9.2 of the summary by the Particle Data Group [jl^]. I have highlighted 
the experimental datapoint from r-decay at ~ 1.8GeV in solid and also included with a solid 
point the value of a s (2.5GeV) ~ 0.2 used in the determination of the gluon condensate by[12|. 



Also shown are two bands derived from estimates of the gluon condensate. The upper band 
corresponds to Narison's[14|, whereas the lower band is obtained from Shifman, Vainshtein 
and Zakharov's[12| value for the gluon condensate. Using ( |7. 1[ ) these estimates imply A^L = 

580 ± 20MeV, respectively A^L = 500 ± 40MeV. To appreciate the quality of the fit note 
that data above Q ~ 20GeV, for which power corrections are negligible, are compatible with a 

(5) 

somewhat larger value of Ai-^. 

to MS 

A recent analysis by Narison using double ratios of moments ||14|| reaches the conclusion 



that 

7 



aSF Hu F Hu) ~ (375 - 400MeV) 4 with a s (1.3GeV) = 0.64±g;?| (7.3) 



<> I Narison 

The advantage of this method compared to the previous one is that the double ratios of 
moments do not depend on the charm mass m c in leading order. The estimate (|7.3|) for 
the gluon condensate is about twice that of|Q. The strong coupling constant a s in this 
analysis was obtained at the optimal scale for the Xc{P\) ~ Xc(Pi) mass-splitting. If a 
perturbative evolution from fi ~ 1.3GeV to fj, ~ Mz is justified, the value for the coupling 
in (|73) corresponds to@ a s (M z ) ~ 0.127 ±0.011 (see also Fig.4). 



With Q the SVZ-value Q for the gluon condensate gives ~ 500 ± 40MeV 



MS 



Using the three-loop evolution [T^] of the corresponding coupling constant (for 



rif = n 
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3) one obtains a s (2.5GeV) ~ 0.325 ± 0.015. Although compatible with other recent 
determinations of a s (2.5GeV) (see Fig. 4), this value is not consistent with a s (2.5GeV) ~ 
0.2 that was used for the sum-rule estimate of the gluon condensate. 



Ag) ~ 580±20MeV. A gain using the three-loop evolution |15| this implies a s (1.3GeV) 



The more recent determination (|7.3|) for the gluon condensate with (|7.1|) corresponds to 



MS 



0.70 ± 0.04, or equivalently, a s (2.5GeV) = 0.356 ± 0.008. While consistent with the 
coupling that was used|14]] to extract the gluon condensate, this value for a s (2.5GeV) is 



somewhat large compared to other determinations of a s (see Fig. 4). Perhaps one should 
note in this context that the second term of the /3-function is almost 40% of the first 
at a s = 0.7. Radiative corrections to the gluon condensate itself (which after all give 
P(g)/2g as the coefficient of (F 2 )) and higher loop corrections on the perturbative side 
of the sum-rules do not appear to be entirely negligible at these large values for a s . The 
quoted error of only 15%[|14| for the gluon condensate is therefore perhaps not a very 



conservative estimate of these systematic errors. 



Using ([74]) and the three-loop evolution |15[ for a s , one can also proceed in the 



opposite direction and estimate ^-j^s an< ^ ^ us ^ ae § mon condensate from a s at scales 
fi < 2.5GeV. This determination of the gluon condensate is only as accurate as the 
estimate for the asymptotic scale, which suffers from the logarithmic relation between the 
latter and a s . Probably one of the best experimental determinations of a s at relatively 
low energies is the one from r-decay. The central value of a s (m T ) = 0.370±0.033 includes 
an estimate of hadronic power corrections to the decay. Neglecting these would, however, 
change a s {m T ) by less than the quoted error |T3]. 



a s (m T ) = 0.370 ± 0.033 corresponds to = 450 ± 80MeV and with Q) implies a 



gluon condensate of {^-F 2 j ~ 300 ± 60MeV. Although only a very rough value for the 
gluon condensate, it agrees within errors with the one obtained by Shifman, Vainstein and 
Zacharov (but for a considerably larger coupling). It is also compatible with Narison's 
lowest estimate. As Fig. 4 shows, this state of affairs is not conclusive, but I consider 
it encouraging in view of the fact that the errors on a s and (^^F 2 ^ are almost wholly 
systematic. 



In summary: the exact relation ( |6.6|) between the asymptotic scale of the modified 
minimal subtraction scheme and (9^) for an SU(3) gauge theory with three quark flavors 
is in reasonable agreement with current estimates of the gluon-condensate and a s . 

Note finally that A-j-f- ~ 500MeV (compatible with Shifman, Vainstein and Zakharov's 
value for the gluon condensate and a s from r-decay) in ( |6.7| ) corresponds to a vacuum 
energy density of about — 455MeV/fm 3 . Assuming that the 5'?7(3)-symmetry is essentially 
restored within a sphere of radius ~ 0.8fm inside a nucleon, one naively obtains M ~ 
970MeV for the mass of such a "bubble". 
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8 Appendix A: Calculation of the two-loop contribu- 
tion to the effective action 



In the gauge 5 = 1, p = the only two-loop contribution to the effective action is the one 
shown graphically in Fig. 2. For a sufficiently large symmetric D = 4 — 2e dimensional 
torus Fig. 2 is the diagrammatic representation for the expression 

(lK(p; j)D ed (q; 7) - D c *{p- 0)D £d (q; 0)) (8.1) 

where g = gp~ e is the dimensionless coupling constant and L is the linear dimension of 
the symmetric torus. In ( |8.1| ) D^ b u (k) is the tree level gluon propagator 

DfM = 5 a \5^k 2 + (a- l)k^k u )/k 4 (8.2) 

and D ab (k; 7) denotes the tree-level ghost- ant ighost correlator for fixed moduli 7. Defining 
the matrix (7)^ = j ahc ^ c and using the integral parametrization 

D ab (k^) = (5 ab k 2 + r^T 1 = rd\e- Xk2 (e-^) , (8.3) 

JO \ / ab 

for the ghost propagators, the D-dimensional momentum integrals in ( |S.1|) can be sepa- 
rated from the color-summations. The expression for Y 2 ~ loop becomes, 

r 2 -'°° p (7) = y A^) D J o °° d\t | o °° dX 2 |/ dac (e~ X ^) cf f fae (e~ A2 % - f dac f cad } 

x / ^yl^ )4 ((P - tfPV + (« - !)(P 2 - M)(PQ ~ I*)) e- Alp2 ~ W (8.4) 
The momentum space integrals can be performed in D > 2 dimensions 



2 



d 2 , <9 2 t a-lfd d \ 2 « + l ^ I f d D pd D q _ Xip 2_ X2q 2_ x(p _ q) 2 



d\ l d\ + d\ 2 d\ + 2 [dXi dX 2 ) 2 «9A 2 ji (2vr/i) 2De 

73 _ 1 _«(£)_ 3) AiA 2 , S5) 



D-2 (Attii 2 ) d (X 1 X 2 ) d / 2 (X 1 + X 2 ] 



2 



I next evaluate the color trace in (|8.4|) in terms of the (real) eigenvalues ej, % = 1, . . . , n 
of the hermitian matrix 27 = it a, y a . Note that 

(e~ xS/ ) ab = -2TrtWe-^ (8.6) 
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with the convention that the anti-hermitian generators t a of the fundamental represen- 
tation of SU(n) are normalized to Tr t a t b = — \8 ab . Repeatedly using the completeness 
relation 

mi = - \ 5 iM) (8.7) 

of the generators in the fundamental representation of su(n), the color trace in ( |8.4|) can 
be rewritten as follows: 

jdac ffae ( e -A a ^ ^ = 4Tr i«]e*iVe- A % [t f , f]e A2 ¥e" A ^ 

= -2Tr [e Xl H f e- Xl \t d }[e X2 H d e- X2 \t f ] 

= Tre A2 ^Tr e Al Ve" (Al+A2) V + Tr e'^Tr e~ Al ¥e (Al+A2 ' V - 2Tr e~ x ^t f Tr e x% V 

— 1 — -Tr e A2 ^Tr e Al ^Tr e ~( Al+A2 ^ — — Tre _A2 ^Tr e~ Al ^Tr e ( ' Al+A2 ^ 
2 2 

= cos(Ai(ej - e k ) + A 2 (e 3 - - e k )) , (8.8) 

where the symbol J2ijk denotes the sum over all triplets (i,j,k) except those with i = 
j = k. The first term in the expansion of the cosine does not depend on the eigenvalues 
and corresponds to the diagram shown in Fig.2c. This 7-independent term of the effective 
action can be absorbed in the normalization M of ( [3.1|) and has been subtracted in ([8.1|) 
since it is of no physical interest. Using (|8.8|) , the color trace in (|8.4j) expressed by the 
eigenvalues of is, 



f dac ( e -A^ //« (e" A2 ^) ed - f dac f cad = £'(1 - cos(A 1 (e l - e k ) + A 2 (e, - e k ))) 

ijk 

= 2^'sin 2 ((A 1 (e i - e k ) + A 2 ( ej - e k ))/2) . (8.9) 

ijk 

When none of the eigenvalues e^, i — 1, . . . , n coincide, the argument of the sine does not 
vanish generically for arbitrary values of Ai and A 2 and regulates the infrared behaviour 
of the integrand in ( |8.4|) . Inserting flS.9|) and (|8.5|) in flS.4|) , the parametric expression for 

r 2-ioop ( -y ~2 ( nD ^-l-«(J-3) w f AidAi [°° X 2 d\ 2 sin 2 (|(Ai^ fc + X 2 v jk )) 
[1) " 9 [fi ] Vte\D - 2) f k Jo A f 12 Jo A ? /2 (Ai + A,)^ 

2 nC D-l-a(I>-3)^/ yi p ^ sin 2 (|A(:rt^ + (1 - x)v jk )) 



~ ° 2{fiL)D 16vr2(D-2) ^Jo dX l dX \ D -\x{l-x))iD-2)/2 

(8.10) 

where V[ m is the dimensionless difference of eigenvalues defined in ( |3.4|) . Feynman's 
parametrization for the integrals shows that the overall degree of divergence of the graph 
is logarithmic, as could have been expected. The integration over A in ( |3.10|) is readily 
performed in D < 4 dimensions and we are left with an integral over a single Feynman 
parameter 

i—w - »)° J ;:,,; ( -V' "(f)^») 
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X 



(xv ik + (1 - x)v jk 



,2\ {D~2)/2 



{x{l — x)) J 

where the symmetry % «-> j under the sum was used. Changing the integration variable to 



y = x/(l — x) and introducing e = (4 — -D)/2, the integral in ( |8.11| ) can also be written 



1/2 



dx 



(xv ik + (1 - x)v j 



Ik) 



> 2 \ (0-2)/2 



1 y 1 £ dy 



v ik + — 

v y , 



io \ x(l — x) J Jo (1 + y) 2 

For e — > 0+ (|8.12| ) diverges as v 2 k /e. To isolate this divergence of (|8.12|) , consider 



;s.i2) 



dy 



l-e 



1 1 



M 2 ^)] + £ o^ + ln ( 4 ) + (Hv%)) ln(4e^)] +. . . 

(8.13) 



o (i + y) 2 V y J ]k V 2 ~~ X ~~ JK ' 1 ' 2 L 6 
and define the for e > —1 and arbitrary real values of a and b regular integral 
■= f Q 1 T T ^y{(ya 2 + b 2 /y + 2aby-c-(by y ) 1 -c} 



dy 



o (1 + y) 



: (ya 2 + 2ab) 



dy 



o (1 + yf 



(ya 2 + 2ab) In \y [a + - H In 



b + ya 



+ 0(e 2 



ab+ (In 2- -)a 2 + 0(e) 



(8.14) 



In gauges a^3 the term of order e in ( 8.14Q gives a finite contribution to r 2 ~' oop in D = 4 
and would have to be computed. For a = 3 - the critical gauge that will ultimately interest 
us - the two-loop effective action simplifies greatly: the factor D — 1 — a(D — 3) in ( |8.11| ) 
in this gauge is of order e. The term of order e in (|8.14| ) thus gives a finite contribution 
to the 4-dimensional effective action proportional to (a — 3). This (finite) contribution 
will be of no further interest and for the sake of brevity will not be given in the form of 
more elementary integrals Q. 



Making use of the symmetric summation over eigenvalues, 
lead to, 

1- r / ' 01 1_e 

L y £ dy / v jk 

o (i +y) [\ y 
l l 



ijk 



ijk 



2n y: 



l<i<j<n 



-\n(v 2 A + 0(e) 



l4|) and ( gl3l) together 

dy T ^%/y) l - £ 

(8.15) 



o (1 + y) 



Expanding ( |8.11| ) for e = (4 — D)/2 — > + using ( |8.15| ), the 2- loop contribution to the 
effective action for a symmetric torus with volume [fiL) D ~ oo gives (|3.9|) of the main 
text. 



4 algebraic programs such as MATHEMATICA readily express the term proportional to e in ( |B.14|) in 
the form of dilogarithms and products of logs. 
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